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Abstract 

In the canonical light-front QCD, the elimination of unphysical gauge degrees of 
freedom leads to a set of boundary integrals which are associated with the light-front 
infrared singularity. We find that a consistent treatment of the boundary integrals leads 
to the cancellation of the light-front linear infrared divergences. For physical states, 
the requirement of finite energy density in the light-front gauge (^4+ = 0) results in 
equations which determine the asymptotic behavior of the transverse (physical) gauge 
degrees of freedom at longitudinal infinity. These asymptotic fields are generated by 
the boundary integrals and they are responsible for the topological winding number. 
They also involve non-local behavior in the transverse direction that leads to non-local 
forces. 



1 Introduction 



Quantum chromodynamics (QCD) was initially proposed as a strong interaction field theory 
in light-front coordinates, motivated by light-front current algebra Q. In recent years, the 
search for nonperturbative solutions of QCD has led to an extensive exploration of light-front 
field theory (LFFT). The main attractions for studying nonperturbative QCD in light-front 
coordinates, called light-front form by Dirac [Q], are that Q: (1) boost invariance in LFFT is 
a kinematical symmetry, which is important in the study of composite systems, particularly 
the hadrons in QCD; (2) LFFT is a relativistic field theory with nonrelativistic structure 
so that the relativistic bound state equations are reduced to Schrodinger-type equations 
from which the nonrelativistic quark model may find its justification in QCD in light-front 
form; (3) the positivity of the longitudinal momentum {k + > 0) in light-front Hamiltonian 
field theory implies that the light-front vacuum consists only of particles with longitudinal 
momentum k + = 0, which may simplify the QCD vacuum structure. These properties 
provide a hope to solve QCD in light-front form for hadrons. 

The earliest systematic formulation of the light-front QCD (QCD in light-front form 
with light-front gauge, referred to simply as LFQCD hereafter) was given about sixteen 
years ago [||. For applications of perturbative LFQCD see ref.||. In order to understand 
basic nonperturbative relativistic bound state problems, in the last few years many works 
on LFFT have mainly focused on various 1 + 1 field theory models, and some on the 3+1 
Yukawa model and QED pi . One main obstacle in extending the study to nonperturbative 
LFQCD is that a formalism to address simultaneously the major difficulties of QCD in light- 
front form is still not in place. These difficulties include the renormalization problem (even 
in perturbation theory), the confinement problem, and the problem of the QCD vacuum 
and dynamical chiral symmetry breaking. 

The renormalization problem in the study of relativistic bound states in LFFT has 
several aspects. Since power counting is different on the light-front [||, there are additional 
ultraviolet divergences in LFFT, compared to the instant form Q. The additional ultraviolet 
divergences have received some attention recently in the context of the relativistic bound 
state problem in the 3+1 light-front Yukawa model ||. LFQCD also contains severe light- 
front infrared divergences. The resolution of the light-front infrared divergence problems is 
not complete even in perturbative LFQCD. Issues arising from the possible mixing of the 
ultraviolet and infrared divergences in the relativistic bound state problems in LFFT have 
not been addressed so far. 

Understanding confinement is crucial for building hadronic bound states in QCD. In the 
present canonical LFQCD, the associated Hamiltonian contains a linear potential between 
color charges only in the longitudinal direction, which does not provide a confinement 
mechanism for quarks and gluons in 3 + 1 dimensions. Therefore, it may not be suitable 
for describing low-energy hadronic structure. Based on light-front power counting, Wilson 
recently proposed a formalism to construct a confining light-front quark-gluon Hamiltonian 
for LFQCD ||. Wilson suggested that a starting point for analyzing the full QCD with 
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confinement in light-front form is the linear infrared divergence (i.e., \/k +2 singularity in 
momentum space). The counterterms for the linear divergence, which may be constructed 
from light-front power counting rules, can involve the color charge densities and involve 
unknown non-local behavior in transverse direction. It is tempting to identify these terms 
as the source of transverse confinement. However, the analysis is not yet complete and a 
scheme for practical calculation has yet to be developed. 

Dynamical chiral symmetry breaking is another important issue in the study of QCD 
for hadrons. In instant form, dynamical chiral symmetry breaking is associated with a 
nontrivial vacuum through the Goldstone mechanism. In LFQCD, the vacuum is trivial 
when the k + = sector is ignored. Therefore, it seems to be natural to argue that in order 
to obtain a nontrivial vacuum, one has to solve the k + = modes Solving the k + = 
modes may provide us with mechanism for spontaneous chiral symmetry breaking. Yet, 
the k + = sector is singular and is very ambiguous. This singularity may exist even in 
free field theory. Thus, it is not clear whether the nontrivial structure of LFQCD must be 
associated with the k + = modes. Furthermore, by involving the k + = sector, the main 
advantage of LFQCD that simplifies nonperturbative bound states is lost, and therefore 
there is no strong reason why we should study nonperturbative QCD in light-front form. In 
fact, dynamical symmetry breaking can be manifested in different ways in different frames. 
It may be more attractive if we could formulate LFQCD with a trivial vacuum such that 
the dynamical breaking of chiral symmetry is manifested explicitly via effective interactions. 
However, the attempt in this direction has not yet started. 

All these problems mentioned above are essential and should involve non-abelian gauge 
degrees of freedom in QCD. We are still unable to solve QCD at the moment. As a start- 
ing point, we shall address in this paper the problems of light-front linear infrared diver- 
gences and the associated nontrivial aspects, based on a canonical quantization approach 
to LFQCD. We hope that these discussions will provide some insight for solving QCD in 
light-front form in the future. 

We apply the conventional canonical procedure[10| to QCD in light-front form. It turns 
out naturally that QCD is a generalized Hamiltonian system where the first-class gauge 
and quark constraints emerge explicitly in the Lagrangian. As is known, in the light-front 
gauge these first-class constraints become solvable first-order differential equations, and 
are used to eliminate unphysical degrees of freedom to all orders of the coupling constant. 
However, the gauge constraint equations contains a set of boundary integrals at longitudinal 
infinity for the longitudinal color electric fields [see eq.(15) and the following discussions]. 
These longitudinal boundary integrals are the color charge density integrated over the longi- 
tudinal space (x~) and are associated with the light-front infrared singularity. The resulting 
LFQCD Hamiltonian contains a boundary term proportional to these boundary integrals 
which are overlooked in previous investigations of light-front gauge theory. 

We find that in perturbation theory the boundary integrals serve to remove linear in- 
frared divergences in loop integrals. Removing the linear infrared divergences in LFQCD 
is a serious problem that has not been solved completely. In usual Feynman theory of per- 
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turbative LFQCD, by use of the gauge fixing term, one can derive the gauge propagator 
involving l/k + singularity. Beyond the leading order calculation, this singularity leads to 
linear infrared divergence in the principal value prescription. In x + -ordering Hamiltonian 
perturbation theory, the linear infrared divergences emerge even in tree-level and one-loop 
diagrams. By including the boundary term in the Hamiltonian, we obtain a consistent 
distribution function for the product of two principal value prescriptions, from which the 
linear divergences in loop integrals are removed by the same divergences in the instanta- 
neous interaction. This finding is useful for perturbative LFQCD calculations in high-energy 
processes. 

The relevant boundary integrals in the Hamiltonian formulation of axial gauge were 
indeed pointed out first by Schwinger in 1962 [12]. Due to the different structure between 
LFFT and the field theory in instant form, the consequences from the boundary integrals 
we study in this paper have not yet been realized explicitly in axial gauge. One of the 
important differences is related to the QCD vacuum. In instant form with axial gauge, the 
QCD vacuum cannot be simple. In LFQCD, generally the vacuum should also be nontrivial 
because of the k + = modes. However, the choice of antisymmetric boundary conditions 
for field variables at longitudinal infinity excludes the k + = modes. In this case, the 
LFQCD vacuum remains trivial as the bare vacuum, and thus the nontrivial QCD structure 
must be carried purely by the boundary behavior of gauge fields. 

For physical states, the requirement of finite energy density results in asymptotic equa- 
tions for transverse (physical) gauge fields at longitudinal infinity. These asymptotic gauge 
fields are generated by the boundary integrals. It is these asymptotic fields that determine 
the gauge field configurations for the non-vanishing winding number associated with the 
topological solutions of QCD. Therefore, although the LFQCD vacuum is trivial with our 
choice of the boundary conditions, the nontrivial behavior of gauge theory is manifested in 
field operators. To the best of our knowledge, this is the first time the nontrivial behavior 
of gauge fields is explored in light-front form with light-front gauge A£ = 0. It is also the 
first attempt to address the nontrivial structure with trivial vacuum in QCD, which seems 
to be possible only in light-front form. 

Furthermore, the asymptotic gauge fields at longitudinal infinity which are generated 
by the boundary integrals not only involve the color charge densities in transverse space 
but also involve non-local behavior in the transverse direction. We find that by replacing 
the nontrivial boundary condition with a trivial one for the transverse gluon fields, many 
transverse non-local potentials are induced by the boundary integrals. These potentials are 
responsible for the nontrivial QCD behavior and therefore may lead to quark and gluon 
confinement. This possibility will be explored in further investigations. 

The paper is organized as follows. In section 2, a canonical procedure for LFQCD is 
studied where we focus on the problem of the boundary conditions in solving the light- 
front constraints. In section 3, the roles of boundary integrals are explored in detail. Some 
remarks are made for relevant problems in section 4. Finally, the paper includes two appen- 
dices. In appendix A, we discuss canonical quantization of LFQCD by use of the rigorous 
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phase space structure 13] rather than the Dirac procedure [14|. In appendix B, we demon- 
strate the cancellation of linear infrared divergences at the tree level in qq scattering and 
in one-loop diagrams for the quark mass correction. 

2 Canonical formulation and boundary condition 

We start from the QCD Lagrangian 

£ = -^Tr(F^F^) + - m)^, (1) 

where = d fl A u -d u A^-ig[A' 1 , A% A» = J2 a A a T<1 is a 3x3 gluon field color matrix and 
the T a are the generators of the SU(3) color group: [T a , T b ] = if abc T c and Ti(T a T b ) = \5 ah . 

The field variable tp describes quarks with three colors and Nf flavors, = \ —igA 11 
is the symmetric covariant derivative, and m is an Nf x Nf diagonal quark mass matrix. 
The Lagrange equations of motion are 

Vr + 9f abc A b ^ + g^T a ^ = 0, (2) 
(i^d" - m + 91„A^ = 0. (3) 

The light-front coordinates are defined as: = x° db x 3 , x l j_ = x i (i = 1,2), where 
x + is chosen as the "time" direction along which the states are evolved, and x~ and xj_ 
become naturally the longitudinal and transverse coordinates. The inner product of any 
two four-vectors is then a^W- = -(a+ft~ -|_ a _ ft+) — a± -b±, and the time and space derivatives 
(cP = ^f-) and the 4-dimensional volume element are given by d + = 2^^, d~ = 2t^, 
d % = — -M^, and d 4 x = ^dx + dx^d 2 x±, respectively. 

Naively, the canonical theory of QCD in light-front form is constructed by defining the 
conjugate momenta of field variables {A£(x), ifi(x), i/j(x)} as 

dC 1 i 

^t(x) = g ( g_^ t) = -'j7 7^ = -^+(x), (6) 

where the fermion spinor in light-front form is divided into ip = ^+ + tp-, V"± = A-±tp with 
A± = ^7°7 ± . Following a similar procedure in instant form described by Faddeev and 
Slavnov[10| for gauge theory, we separate the time derivative terms from the Lagrangian, 



-|A-C a + i(^C + cV_)}, (7) 
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where 



W = 2 + B- 2 ) + I • + sAl) + 0m}i/>- + /i.e.} 

+ (8) 



and 

Ca = \{d + E- + gf ahc AtE-) - (d'Ei + gf^J&ft + g^ + T^ + (9) 
C = (id + + gA + )ip_ - (ia± ■ d± + ga ± ■ A± + j3m)^ + . (10) 

In eq.(8), we have defined B~ = F^ 2 as the longitudinal component of the light-front color 
magnetic field. 

The reason for writing the Lagrangian in the above form is to make the Hamiltonian 
density and also the dynamical variables and constraints manifest. In eq.(7), the first term 
contains all the light-front time derivative terms. From the definition of eqs.(4-6), it immedi- 
ately follows that only the transverse gauge fields A\ and the up-component quark fields '(/'+ 
and are dynamical variables. The second term in eq.(7), 7i, is a Hamiltonian density. It 
contains three parts, the first part involves the light-front color electric and magnetic fields; 
the second, the usual quark Hamiltonian with coupling to the gauge field, and the last a 
surface term. Besides the kinetic term and the Hamiltonian density, eq.(7) also contains 
an additional term. This is a constraint term which indicates that the longitudinal gauge 
field A~ and the down-component quark fields if)- (ipL.) are only the Lagrange multipliers 
for the constraints C a ,C (C') = 0. These constraints arise from the definition of canonical 
momenta in the light-front coordinates and are consistent with the Lagrangian equations 
of motion. The gauge field constraint, C a = 0, is in fact the light-front Gauss law which 
is an intrinsic property of gauge theory. The fermion constraint, C (C') = is purely a 
consequence of using the light-front form. 

The existence of constraint terms simply implies that QCD in the light-front form is 
a generalized Hamiltonian system[ll]. These constraints are all secondary, first-class con- 



straints[| in the Dirac procedure of quantization [14]. To obtain a canonical formulation 
of LFQCD for non-perturbative calculations, we need to explicitly solve the constraints, 
namely to determine the Lagrange multipliers, to all orders of the coupling constant. Gen- 
erally, it is very difficult to analytically determine the Lagrange multipliers from the con- 



straints C a ,C = since they are coupled by A£. Only in the light-front gauge [15, p" 



A+(x) = A° a (x)+Al(x) = 0, (11) 



1 In light-front field theory, there always exist the so-called primary, second-class constraints due to the 
fact that the light-front Lagrangian is linear in the first-order i + -derivative and therefore the canonical 
momenta are functions of field variables. These constraints are not real constraints in the generalized phase 
space quantization procedure and are easily handled. We will put these discussions in the Appendix A. 
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are these two constraints reduced to solvable one-dimensional differential equations: 
\d + E~ = &EI + g(f abc AlEl - ^ + T a ^ + ) = G a 

(12) 

id + ijj- = (ia± • d± + ga± • A± + (3m)tp + 
In order to solve eq.(12), we have to define the operator l/<9 + . In general, 

(^■^j f(x~,x + ,x ± ) = J dx^e(x~ - Xi)f(xi,x + ,x ± ) + C(x + ,x ± ) (13) 

where e{x) = 1,0,-1 for x > 0,= 0, < 0, respectively, and C(x + ,x±) is a x~ independent 
constant. However, since the canonical conjugate of transverse gauge field in LFQCD is a 
dependent variable [E l a = —^d + A l a , see eq.(4) with A+ = 0], one has to impose a priori a 
boundary condition for A\ in order to derive the canonical commutation relations for the 
physical field variables. It has been shown |l6|, [lTj that the suitable definition of \/d + which 
uniquely determines the initial value problem at x + = for independent field variables is 
C(x + ,x±) = 0. This corresponds to choosing an antisymmetric boundary condition for 
field variables in the longitudinal direction. 

Using eq.(13), we can explicitly express E~ in terms of transverse gauge fields A\ and 
the independent light-front quark field ip+ from the gauge constraint in eq.(12), 

E-{x) + d i A i a {x) = - 9 ~ r dx / -e(x--x'-)(f abc Ald + Al + 24;iT a 4> + ) + C a (x + ,x ± ), (14) 

4 J— oo 

here, E l a = — ^d + A l a has been used. To uniquely determine their initial values at x + = 0, 
we require that the E~ and A l a satisfy antisymmetric boundary conditions at longitudinal 
infinity, namely, C a {x + ,x±) = 0. As a result, 

E-(x) = - 9 - r dx'~e(x- - x'-)(f abc A\d + Al + 2^T> + ) 

4 J —oo 

X ' = =°° -9Va-=±oo T 7 r dx-(f abc Aid+Ai + 2^T> + ). (15) 

4 J —oo 

Since E~ satisfies now an antisymmetric boundary condition, its boundary values at 
longitudinal infinity are completely determined by the second equality in eq.(15), where the 
second term is boundary integrals over x~ for the color charge densities. These integrals 
are the source of light-front infrared singularity. We call them the longitudinal boundary 
integrals, or simply the boundary integrals. 

By using the identity) 12], 

1 f x / 2 _____ 1 

- / dx e(x — x' )e(x — x" ) = — \x — x" I H — A, (16) 

2 J-x/2 2 
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where the parameter A denotes the distance between two boundary points in the longitudinal 
direction, the color electric field energy in the Hamiltonian becomes 



1 roo 1 roo f 

H E = - dx-d 2 x ± (E-) 2 = - / dx-d 2 x L (d'Al) 2 

^ J —oo ^ J —oo I 

1 r°° r 

+ - / dx'~ gd l A l a e{x~ - x'"){f abc A l b d + A l c + 2^\T a ^ + ) 

2 J —oo 

fl 2 r-OO 

- 9 - / dx'-(f abc Ald + Ai + 2^T> + ) 

o J-oo 

\x~ - x'-\{f ade A 3 d d + A{ + 2^T>+) | 

+ (hm A) ^ J d 2 x ± y^dx-(r hc Ald + Al + 2ip{T a i> + )}\l7) 

In the above equation, the last term (a boundary term) involves the boundary integrals 
and is associated with the infrared divergence in the light-front instantaneous interactions. 
As we will discuss later, in perturbation theory, this term is regularized by the distribution 
function of the product of two principal value prescriptions and leads to the cancellation 
of the light-front linear infrared divergences. For physical states, the requirement of finite 
energy density results in the asymptotic equations for the transverse gauge fields which show 
that the asymptotic transverse gauge fields do not vanish at longitudinal infinity and are 
generated by the boundary integrals. The non-vanishing asymptotic transverse gauge fields 
determine the nontrivial QCD structure in light-front form. Thus, the boundary integrals 
can inherently affect QCD dynamics. 

The Lagrange multipliers in eq.(7) can be determined easily now. The Lagrange multi- 
plier is the solution of the quark constraint in eq.(12), 

ip-(x) = — — / dx'~d 2 x'_ L e(x~ - x'~)5 2 (xj_ - x' ± ) 

4 J —oo 

x [ a± ■ (id' ± + gA ± (x')) + 0m]^+{x'). (18) 
The Lagrange multiplier A~ is obtained from the definition E~ = — ^d + A~ and eq.(15), 

1 f°° 

= -o / dx'-e(x- -x'-)E-(x + ,x'-,x ± ) (19) 

^ J —oo 

For this solution, the first surface term in eq.(8) vanishes. Moreover, it is reasonable to 
assume that the transverse color electric fields E l a as well as A % a vanish as 0(r~ 2 ) and 
0(r _1 ) at r = \x±\ — > oo because the gauge freedom is totally fixed at the transverse 
infinity. Hence the other surface term in eq.(8) vanishes as well. 

After the determination of the Lagrange multipliers, the LFQCD Hamiltonian is given 
simply by 

H = J dx~d 2 x ± S^(E- 2 + B- 2 ) + Vj>i(^i + gA*) + 
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/ 



dx-d 2 x± + gf^AiAi&Ai + y^r^AiAiAiAi 




-z-0+{ai(i9i + gA { ) + f3m}e(x~ - x'~){a? L {id> ± + gA j ) + /?m}V>. 




/"OO 

/ dx'- {f abc Aid + Ai + 2^r> + ) 



a;" - x'-\{f ade A j d d + A{ + 2^^ a V+)} 




(20) 



A detailed procedure to quantize the above formulation is presented in appendix A. This 
Hamiltonian contains the naive canonical LFQCD Hamiltonian plus a boundary term which 
is the square of the boundary integrals, as a result of eq.(15). Choosing antisymmetric 
boundary conditions in the longitudinal direction in LFQCD has the following advantages: 

1) . With any other boundary condition, eq.(14) contains an arbitrary ^"-independent 
function. Such an arbitrary function leads to ambiguities in formulating LFQCD. Only with 
an antisymmetric boundary conditions, is this arbitrary term zero and formally LFQCD can 
be completely defined. 

2) . For the definition eq.(13) with C(x + ,x±) = 0, all field variables in LFQCD satisfy 
antisymmetric boundary conditions at longitudinal infinity except "0+( x )> whose boundary 
condition is not specified. However, the equation of motion for if)+(x) contains 

[see Appendix A] which forces it to satisfy the antisymmetric boundary condition. As a 
result, the k + = modes are completely excluded in the momentum expansion of field 
variables. Since the LFFT vacuum is occupied only by the k + = particles, with the help 
of antisymmetric boundary condition the LFQCD vacuum is ensured to be trivial as the 
bare vacuum. The important consequence is, as we shall discuss in the next section, that 
in this case the nontrivial structure of QCD is carried purely by field operators. 

3) . By the choice of antisymmetric boundary conditions, the residual gauge freedom in 
A+ = is completely fixed |l8| (also see the discussion later). 

It is worth pointing out that there is a disadvantage for using eq.(13) in perturbative 
light-front Feynman loop integrals. In perturbative theory, eq.(13) with C = leads to the 
principal value prescription for the light-front longitudinal infrared singularity. Although 
this prescription removes severe linear infrared divergences as we will show next, severe 
logarithmic infrared divergences are still present. These logarithmic infrared singularities 
correspond to the "spurious" poles in Feynman integrals, which prohibit any continuation 
to euclidean space (Wick rotation) and hence the use of standard power counting arguments 
for Feynman loop integrals This causes difficulties in addressing renormalization of 

QCD in Feynman perturbation theory with light-front gauge. In the last decade there are 
many investigations in attempting to solve this problem. One excellent solution is given by 
Mandelstam and Leibbrandt, i.e., Mandelstam-Leibbrandt (ML) prescription |20|j , which 
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allow continuation to euclidean space and hence power counting. This prescription has also 
been derived in the equal-time quantization with light-front gauge by Bassetto et al. pl| . 
It has been also shown that with ML prescription, the multiplicative renormalization in the 
two-component LFQCD Feynman formulation is restored j22|. 

In the present paper, we study QCD in light-front equal-x + quantization. Unfortunately, 
ML prescription cannot be applied to equal-x + quantization because ML prescription is 
defined on the boundary condition which involves x + itself [^3| and are not allowed in equal- 
x + canonical theory. Yet, as is pointed out recently by Wilson ||, light-front power counting 
differs completely from the power counting in equal-time quantization. Furthermore, the 
current attempts to understand nonperturbative QCD in light-front form is based on the 
x + -ordered (old-fashioned) diagrams in which no Feynman integral is involved || Thus 
the problem with the power counting criterion in Feynman loop integrals does not affect 
our discussions in this paper. The renormalization of light-front Hamiltonian is an entirely 
new subject where investigations are still in their preliminary stage |9|. 

The main aim of this paper is to show that the boundary integrals play a important role 
in understanding the nontrivial features of LFQCD. The logarithmic infrared divergences 
are completely cancelled in the complete loop diagrams of dynamical processes, as was 
previously shown in the calculation of QCD correction to the scale evolution of hadronic 
structure function up to two- loop |24|. A simple example of such a cancellation in x + - 
ordered perturbative theory is also given for quark mass renormalization in Appendix B-2. 
In our forthcoming papers |25| we will present a detailed discussion on x + -ordered perturba- 
tive loop calculations and light-front renormalization in LFQCD Hamiltonian theory, where 
the logarithmic infrared divergences are again completely cancelled in coupling constant 
renormalization. However, based on light-front power counting, the linear infrared diver- 
gences only involve color charge density and involve non-local behavior in the transverse 
direction, which may be the source of transverse confinement in LFQCD ||. The severe 
linear infrared divergences in LFQCD has not been explored in light-front Hamiltonian. In 
the present paper, we shall focus on linear infrared singularity associated with the boundary 
integrals, which may relate to the nonperturbative aspects of LFQCD in physical states, as 
we will see below. 



3 Role of boundary integrals 

1. Removing linear infrared divergences. In the past decade, applications of LFQCD are 
mostly restricted in perturbation theory. Naively, the boundary term in eq.(20) is ignored 
so that the light-front instantaneous interactions are thought to be linear potentials |l5|, |26fl . 
However, this negligence leads to severe infrared singularities in the perturbation theory. 
To see this clearly, we consider the formulation in momentum space. For the prescription 
of l/<9 + expressed in terms of the integral of eq.(13), the standard Fourier transform leads 
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to the principal value prescription in momentum space as follows, 



(^p) /(O = \ j^^"^ 1 ' " "O/OO 

(^p) /OO = i2 /°° dx'-dx"-e{x- - x'-)s(x'- - x"-)f(x"-) 

f(k + ) = ( 22 ) 

Eq.(22) defines the product of two principal value prescriptions of eq.(21) in terms of the 
distribution function. In this derivation, it follows (see eq.(16)) that the boundary term in 
eq.(20) have been regularized. It is known that eq.(22) leads to linear infrared divergences 
in loop integrals. In order to avoid this divergence, one naively introduces the following 
prescription || 

Vq+) f( x ) = ^2 J_ dx \ x ~ x \f( x ' ) 

-Kw^ + w L w) }(k+ ^w\ fik+l (23) 

This corresponds to the case that the longitudinal boundary term in eq.(16) is ignored. 
Equivalently, the last term in eq.(20) is dropped. Apparently, this prescription removes the 
linear infrared divergence originated from the instantaneous interactions. Unfortunately in 
such a prescription, beyond leading order calculations in Feynman perturbation theory or 
even in leading order calculation in the old-fashion Hamiltonian perturbation theory, the 
product of two principal value prescriptions appearing from three-point vertex either is not 
defined or leads to linear infrared divergences. We shall show that it is the prescription 
of eq.(22) which serves for the cancellation of linear infrared divergences originated from 
the three-point vertex and from the instantaneous interactions. Here we only discuss the 
x + -ordered (old-fashion) perturbative calculations. 

First at the tree level, for example for qq scattering (see appendix B-l), only the linear 
potential leads to a 1/e 2 divergence as k + — > 0. The scattering involving one-gluon-exchange 
is finite due to the principal value prescription. Thus in the naive prescription (23), even 
the lowest order qq scattering amplitude is 1/e 2 divergent. By including the boundary term, 
this divergence is cancelled. In loop calculations, for example for the one-loop correction to 
the quark self-energy (see appendix B-2), the one gluon exchange diagram (which contains 
an integral of l/[£; + ] 2 ) leads to a 1/e divergence; the linear potential is, however, infrared 
finite in the relevant integral of l/[fc +2 ]. Hence, in the naive prescription, loop calculations 
also contain the severe 1/e infrared divergence. In prescription (22), the instantaneous 
interactions are the linear potentials accompanied by the boundary term, which produce 
a 1/e divergence from the integral of l/[fe + ] 2 that cancels precisely the same divergence in 
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the one gluon exchange diagram. Furthermore, the cancellation in the one-loop correction 
of quark-gluon vertex has also been verified |25]| . 

The reason that the linear infrared divergences are removed by using prescription (22) 
can be understood as follows. From eq.(15), the k + singularity originated from the boundary 
integrals. The color electric energy in LFQCD Hamiltonian contains two sources for the k + 
singularity. One is the explicit boundary term, the last term in eq.(20), which is l/k +2 - 
singular. The other belongs to the gluon emission vertex. The resultant gluon emission 
vertex is the first term in the square bracket in eq.(20), which is l/A: + -singular. Therefore, 
in one gluon exchange diagrams, it produces a l//c +2 -singularity, namely the product of two 
principal value prescriptions for the definition of eq.(21). The associated linear divergence 
in loop integrals is the same as that from the l/A; +2 -singularity of the boundary terms in the 
prescription of eq.(22), with a different sign from an energy denominator, and therefore the 
linear infrared divergence is cancelled. Note that in eq.(20) there is another l//c + -singularity 
[in the second term in the square bracket], which comes from the quark constraint [see 
eq.(18)]. Yet, in one-gluon exchange diagrams, it leads to a form l/(pfp2) (pf = p\ + k + ) 
which does not generate infrared divergences. Thus, all linear infrared divergences originate 
from the same source, the boundary integrals. Any negligence of boundary term in the 
Hamiltonian through eq.(23) will lead to unwanted infrared divergences. 

However, the cancellation of the linear infrared divergences in higher order loop-integrals 
(beyond the one-loop diagrams) may also depend on the regularization of ultraviolet diver- 
gences. The cancellation beyond leading order should be true for gauge invariant regular- 
ization. For gauge variant regularization, such as transverse dimensional regularization Q, 
boost invariant cutoff regularization [27] and the explicit cutoff regularization [28] used in 
the x + -ordered perturbative LFQCD, we need to introduce gluon mass counterterms. These 
counterterms break gauge invariance and thereby may also spoil the cancellation of linear 
infrared divergences in higher order diagrams. However, if we set the quark mass m = 
in perturbative LFQCD, the transverse dimensional regularization results in a zero gluon 
mass correction. In this case the cancellation is still satisfied in two-loop diagrams. In deep 
inelastic scattering, one often sets m = in calculating high-order corrections to the scale 
evolution of hadronic structure functions [^| . A more detailed discussion on perturbative 
LFQCD will be presented in a separate paper [25]. For low-energy dynamics, the light quark 
mass is crucial and perturbation theory is no longer useful. Removes infrared divergences 
needs to be treated in an alternative way, which we shall discuss later. 

We may point out that in 1 + 1 LFQCD p9fl , the boundary integral is the color charge 
operator. The corresponding boundary term occurring in the Hamiltonian [30] is then 
proportional to the square of color charge. It is indeed this term resulting in an infinite 
quark mass which is regarded as evidence of quark confinement in 1 + 1 QCD. Explicitly, 
the linear potential does not provide an infinite mass for the quark, as shown above (also see 
ref.(3l|]), but the boundary term adds a l/k +2 singularity (the 1/e divergence) to the quark 
propagator. Since there are no transverse gluons in 1 + 1 QCD to cancel this divergence, the 
boundary term recovers 't Hooft's solution of the infinite quark mass pole|29|]. In physical 
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(zero color charge) states, the boundary term does not contribute to physical observables 
since it is the square of the color charge operators. Quark confinement in gauge-invariant 
states arises purely from the linear potential. This implies that ignoring the boundary 
integral in 1 + 1 QCD may not affect any observable. In 3 + 1 QCD, the existence of 
transverse gluons changes these consequences. 

2. Winding number. The winding number is a topological quantity associated with 
nontrivial structure of U(l) problem in QCD (in Euclidean space, this is the topological 
charge or Pontryagin index associated with the instanton solution of non-abelian gauge 
theory) |3^, 33 1. The second important consequence from the boundary integrals is that 
they determine a non-vanishing winding number. 

As it is known y4+ = cannot completely fix the gauge degrees of freedom. There 
are residual gauge transformations under which the theory is invariant. The generators of 
residual gauge transformations are 

Ra = "2 J dx~ {2d i d+Ai + g(f abc Ald + Ai + 2V> + T> + )} (24) 



As we have shown |18|] the residual gauge transformations generated by R a break the anti- 
symmetric boundary condition for A l a at longitudinal infinity and therefore are not allowed 
with our choice of the boundary conditions. However, by using the antisymmetric boundary 
conditions, the first term in eq.(24) can be integrated out explicitly over x~ and we have, 

R' a = T^Ai\ x - =±00 - 9 -T dx-(f abc Ald + Al + 2^rV+). (25) 

£ J—oo 

It has been verified that for this definition, R' a generate the gauge transformations 
preserving the antisymmetric boundary conditions. But these operators do not commute 
with the LFQCD Hamiltonian (20). In other words, there is no longer an additional gauge 
freedom to choose other A l a such that the resulting Hamiltonian remains invariant. There- 
fore, with antisymmetric boundary conditions, the residual gauge freedom in = is 
completely fixed fi~8|| . 

Furthermore, for physical states, finite energy density requires that the longitudinal 
color electric field strength must vanish in the infinity (A similar requirement was used by 
Chodos in axial gauge |@): 

E-\ x - =±oo = (26) 



or explicitly 



0U* L-=±oo = dx-{f abc Ald + Al + 2^T<V + ). (27) 

4 J — oo 



Eq.(27) is consistent with our choice of antisymmetric boundary condition. Moreover, 
this condition explicitly shows that the transverse gauge fields at longitudinal infinity are 
generated by the boundary integrals.^ Clearly, eq.(27) is satisfied only for physical states. 
In perturbation theory, we cannot use this condition because in perturbative QCD, we 

2 It may be worth pointing out that unlike the 1 + 1 QCD, the boundary integrals in 3+1 LFQCD are 



13 



consider not only physical states but also color non-singlet states for which eq.(27) may not 
be satisfied. Therefore the main effect of eq.(27) should be manifested in nonperturbative 
dynamics. Now we shall discuss how eq.(27) determines the non- vanishing topological 
winding number. 

To make the discussion clear and without any loss of generality, we consider the case of 
zero quark mass only|i4|]. It is well-known that the axial current for iVj-flavor quarks has 
an anomalous divergence, 



d ^ = N fh>^( F ^l (28) 



8tt 2 

where the axial current is j§ = ipj^^ip, and the dual field strength is F^ v = ^e^ uap F cr p. 
From eq.(28), the time derivative of the light-front axial charge is given by 

d'Qb = J dx~d 2 x ± d"{i^ + -f b tp + ) 

-,2 



Nf h>J dx~d 2 x±Tr(F liV F ,J ' v ). (29) 



In obtaining eq.(29), we have used the fact that the other three surface terms of axial 
currents at longitudinal and transverse infinity do not contribute to d~Q^. This is clear 
if we note that j\[x) — > at x± — > ±oo and j§(x) = 2ip_(x)^/ 5 ip_(x) which leads to 
\ x - =00 — js\x-=-oo = due to the antisymmetric boundary behavior of ip_ at longitudinal 
infinity (see eq.(19)). 

The anomaly alone does not imply that the right-hand side (r.h.s.) of eq.(29) must 



be nonzero) 32]. The nonzero contribution of the r.h.s of eq.(29) is given by the gauge field 
configurations determined by eq.(27). This can be seen from the winding number in LFQCD 
defined as the net charge between x + = — oo and x + = oo, 



AQ 5 = \ J dx + d~Q 5 = N f £z J d 4 xTr(F^F^). (30) 

The integration on the r.h.s. of the above equation is defined in Minkowski space. By using 
the identity 

Tr(F^F^) = 43^ (31) 

where 

K* = e^Tr U v d a A p + \A u A a A p } , (32) 

color charge densities in the transverse space and not the color charge operators. The color charge operator 
is defined as follows: 

Q a = J dx-d 2 x ± (r bc Ald+Al + 2iPlT a 4< + ) 
= J dx~d 2 x ± (p°(x- ,x ± ) + pl(x~,x±)) 

where p 9 a is the charge density carried by gluon field and p q a that by quark field in light-front form. Therefore, 
we cannot simply drop the boundary integrals for color singlet states. In other words, removing the boundary 
integrals in eq.(20) implies ignoring the non-vanishing asymptotic A r a fields at longitudinal infinity, and 
therefore losses the possibility to address the nontrivial properties in LFQCD. 
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the r.h.s. of eq.(30) is reduced to surface integrals. In the light-front gauge A+ = 0, the 
second term in K + ' 1 is zero. Thus, the surface integrals at transverse infinity vanish because 
if* falls off as r~ 3 for r = \x±\ — > oo. Meanwhile, since x + and x~ are symmetric, it is 
reasonable to use the same boundary condition (i.e., antisymmetric boundary condition) 
for A l a at x + = ±00. Therefore, the contribution from the surface integral at x + -infinity 
vanishes as well. At x~ = ±00, the surface integral contribution for the first term of K~ is 
also zero due to the antisymmetric boundary condition for the A\ and A~ fields. Finally, 
eq.(30) is reduced tcQ 

AQ 5 = -N f ^- / dx + d 2 x ± Tr (A - [A 1 , A 2 ]) 

IT J ' »-=— 00 

= -2iV/^2 / dx + d 2 x ± Tt(A~[A\A 2 ]), at x~ = 00, (33) 

Here we have used again the antisymmetric boundary condition of A\ and A~ at longitudinal 
infinity. Eq.(33) shows that a non- vanishing AQ5 is generated from the asymptotic fields of 
A l a ,A~ at longitudinal infinity, which are generated by the boundary integrals by eq.(27). 
The non-zero A~ | x - = _|_ 00 is determined by transverse gauge fields (from Eq.(19)), 

1 f°° 

A~(x + ,x±)\ x - =±oa = ±- dx-x~G a {x + ,x-,x ± ). (34) 

£ J — 00 

Thus, the boundary integrals are essential for the non-vanishing winding number in QCD. 

The above derivation shows that for an antisymmetric boundary condition, although 
the LFQCD vacuum is trivial, the nontrivial QCD structure is switched to the field opera- 
tors. This structure is manifested in the asymptotic behavior of transverse gauge fields at 
longitudinal infinity and is explicitly associated with the boundary integrals. The trivial 
vacuum with nontrivial field variables in the present formulation of LFQCD may provide a 
practically useful framework for describing hadronic states. We now turn to this discussion. 

3. Non-local potentials in the transverse direction. One of the nonperturbative ap- 
proaches to solve bound states in LFFT is the Tamm-Dancoff approach, which truncates 
the Fock space to be a few-body state space |35]] . Such an approach becomes practically 
applicable only when the vacuum is trivial. We have given a realization of a trivial vacuum 
in this paper. However, to address hadronic bound states, the existence of nontrivial poten- 
tials, namely, confinement potentials, is crucial. An explicit construction of such potentials 
from QCD is still lacking. The LFQCD Hamiltonian contains linear potentials only in the 
longitudinal direction (see eq.(20)). Quark and color confinement certainly requires similar 
potentials in the transverse direction as well. We suggest that these nontrivial potentials in 
LFQCD might hide in the condition of eq.(27). 

From eq.(27) we see that the asymptotic A l a fields at longitudinal infinity are propor- 
tional to the color charge density in transverse space and also that they involve non-local 



3 As we have seen here, this solution is obtained by use of eq.(28). It would be very interesting if one 
could directly derive AQ5 or eq.(29) from the LFQCD Hamiltonian eq.(20). 
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behavior in the transverse direction (induced by the transverse derivative). Intuitively, we 
may separate the transverse gauge potentials into a normal part plus a boundary part, 

Ai = A i aN + Ai B (35) 

where 

^Kin\x~=±oo = > 9 l A l aN ^-=±00 = 0, (36) 

VA\ B \ x - =±00 = t| (P*(*±) + P q a (x±)). (37) 

In eq.(38), p a {x±) denote the color charge densities integrated over x~ , where the definition 
of the color charge densities is given in footnote 2. The conditions of eqs.(36) and (37) do not 
uniquely determine the separation of eq.(35). Generally, there are two types of separation for 
eq.(35). One is to consider A l aB the long-distance fields generated by the boundary integrals 
and A l aN the short-distance fields determined by free theory. Thus, A l aB correspond to the 
gauge field configuration for the non-vanishing winding number discussed in the previous 
subsection. In this case, if we are only interested in the low-energy dynamics, the effect 
of the A l aN fields may be ignored. Thus, it is very attractive but is also very difficult to 
analytically find the A l aB . Another possibility is to choose a simple solution for the A l aB 
that satisfy eq.(37). In this case, the A l aN have the trivial boundary condition eq.(36) but 
are not determined by free theory. The Hamiltonian is then expressed only in terms of the 
A l aN , and the boundary behavior of transverse gauge fields are replaced by the effective 
interactions. A convenient choice for A l aB which satisfy eq.(37) is 

Al B (x) = -^Jdx'~dx' l e(x'-x^e(x l -x' l )(pi(x')+pUx')) , i = 1,2. (38) 

Substituting the separation of eq.(35) with (38) into the LFQCD Hamiltonian, we obtain 
a new Hamiltonian in terms of A l aN that contains many effective interactions induced by 
eq.(27). All these effective interactions involve the color charge densities and involve non- 
local behavior in both the longitudinal and transverse directions. One of the lowest order 
interactions, for example, is given by 

/oo 
dx 1 dx~ dx> dx'~ dx' 3 rf 3 {d % p q a {x~ , x l , x j )\x~ - x'~\ 

u "°° ' 

\x j -x' j \d i pl{x , -,x i ,x ,: ')} (39) 

where rf 3 = 1 (0) for i ^ (=) j. Hence, eq.(37) leads to numerous many-body non-local 
color charge interactions which are functions of boundary integrals, and which may lead to 
confinement. 

As we have mentioned in the introduction, Wilson recently proposed a formalism to con- 
struct a confining light-front quark-gluon Hamiltonian for LFQCD ||. Wilson suggested 
that a starting point for analyzing the full QCD with confinement in light-front form is 
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the linear infrared divergence (i.e., l/fc +2 singularity in momentum space). Our proce- 
dure for constructing effective interactions from the boundary integrals is indeed associated 
with the linear infrared divergences. However, Wilson's approach is totally different from 
what we have discussed here. His analysis based on light-front power counting shows that 
possible confining potentials may be obtained from the counterterms of the linear infrared 
divergence. In QED, the counterterms for infrared divergences are forbidden because the di- 
vergences arise from integration over the square of the electron scattering amplitude rather 
than integrals over the amplitude itself and it is only the latter that can be regulated by 
counterterms. However, in QCD, since quarks and gluons serve only as constituents, there 
are no scattering cross section for them, and there exist counterterms for the linear infrared 
divergence. These counterterms involve the color charge densities integrated over x~ and 
involve non-local behavior in the transverse direction. However, the non-local structure of 
the counterterms is unknown since the power counting itself cannot determine it. They 
also violate longitudinal boost invariance so the coefficients of all the counterterms may 
further be determined by the requirement of boost invariance. Wilson pointed out that it 
is tempting to identify these terms as the source of transverse confinement ||. 

In the present formulation, as we have seen from eq.(27) or (38), the asymptotic gluon 
fields which induce effective interactions are proportional to color charge densities integrated 
over x~ (the boundary integrals) and also involve non-local behavior in the transverse 
direction. Furthermore, the non-local behavior of effective interactions in the transverse 
direction is also determined by eq.(27) or (38). Thus, we can explicitly construct many 
effective interactions from (38). It can be shown that the effective interactions related to 
fermion part are similar in both QED and QCD but are very different for the gauge part. In 
QED, there is no any effective interaction generated by the boundary integrals that involves 
gauge fields. However, in QCD, there are numerous number of effective interactions which 
are coupled to the non-abelian gauge fields. These effective interactions may be responsible 
for quark confinement since they originate from the nontrivial behavior of the non-abelian 
transverse gauge fields. Yet, it is interesting to see that the analyses based on different 
approaches have the same consequence that non-local potentials in the transverse direction 
in LFQCD are related to linear infrared divergences which have not been paid attention in 
the previous investigations. 

Still the Hamiltonian contains, in principle, an infinite number of many-body interac- 
tions generated by the boundary integrals (or obtained from the counterterms of the linear 
infrared divergences). This is a consequence of the boundary integrals in a non-abelian 
gauge theory due to the existence of nonlinear gluon interactions. It is also true in other 
gauge choices, such as Coulomb gauge [36] or axial gauge |T2|. Practically, as the first step, 
we may only keep two-body interaction terms, such as eq.(41), in the new Hamiltonian. 
Because of the trivial vacuum in the present formulation of LFQCD, using such an approx- 
imate LFQCD Hamiltonian, we can apply the light-front Tamm-Dancoff approach to find 
hadronic bound states, where the bound states contain only a few particles, such as one 
quark-antiquark pair, one quark-antiquark pair with one and two gluons. This is certainly 
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one of the most attractive approaches for low-energy QCD. A numerical investigation along 
this consideration is in progress. 



4 Discussions 

In the previous section, we have discussed some primary properties of boundary integrals 
which we think to be important in understanding LFQCD. In the current investigations of 
LFFT, one of the most active topics is the problem of the k + = modes and the A + = 
gauge. The implications of the boundary conditions in determining the nontrivial behavior 
of gauge theory has, however, been overlooked. In this section, we have some remarks to 
make about the relation of the k + = modes, the = gauge and boundary conditions 
at longitudinal infinity. 

In previous LFFT investigations, much attention has been paid to how to construct a 
nontrivial vacuum from the k + = modes. All attempts have focused on 1 + 1 field models 
0. The motivation for these attempts, as we have mentioned in the introduction, is to 
try to understand spontaneous chiral symmetry breaking in LFFT. In instant form, the 
vacuum is, of course, crucial for hadronic structure since we believe that axial charges Q§ 
create pseudoscalar particles (the lowest bound states in strong interaction region) from the 
vacuum. However, the role of light-front axial charges in hadronic structure is totally differ- 
ent. The success of light-front current algebra in describing low-energy hadronic structure 
is based on the properties of light-front Q5 with a trivial vacuum |37j . In this case, Q§ anni- 
hilate the vacuum so that the vacuum in LFFT itself is not essential in understanding chiral 
symmetry. The importance of light-front Q% lies in their matrix elements between hadronic 
states. These matrix elements are proportional to hadronic decay constants involving pseu- 
doscalar mesons and therefore carry the basic information of hadronic structure. In instant 
form, the matrix elements of Q% in hadronic states with zero momentum transfer are zero 



if one does not make use of the infinite momentum limit 3£]. In other words, the instant 
Q\ itself is practically not useful for hadronic structure except for the Nambu-Goldstone 
picture of spontaneous symmetry breaking, where the important ingredient is the axial 
current. These totally opposite properties of axial charges in light-front and instant forms 
implies that to address dynamical breaking of chiral symmetry in LFQCD, one may need 
to understand the relation between light-front axial charge operators and the Hamiltonian 
operator rather than the structure of vacuum in LFFT. The present work is motivated by 
this consideration. 

Second, the canonical quantization of light-front gauge theory is often considered in a 
box with a periodic boundary condition for gauge fields. In this case, it was argued that 
one cannot use = for the k + = sector because it is incompatible with the periodic 



boundary condition in finite volume |3S]. This argument is true but incomplete. The 
A£ = condition does not totally fix the gauge freedom, and the residual gauge fixing is 
responsible for the non-trivial gauge field configurations which, however, are lacking when 
one imposes a periodic boundary condition [see eq.(33)]. Thus, if one prefers to use a 
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periodic boundary conditions to quantize LFQCD in a box, one has to choose other gauges 
p9| . However, for any gauge fixing other than the A+ = 0, eqs.(9) and (10) show that the 
constraint conditions are extremely difficult to solve except for numerical calculations, such 
as lattice gauge calculations. Our treatment in this paper shows that we can address the 
non-trivial QCD structure by use of A+ = gauge and the trivial vacuum if we take into 
account the residual gauge fixing in antisymmetric boundary conditions for field variables 
at longitudinal infinity. In principle for A+ = gauge, other boundary conditions can also 
be used [40] yet the resulting theory is currently intractable due to the existence of the 
k + = modes. 

Finally, we discuss briefly the difference between LFQCD and the canonical formulation 
of QCD in instant form with axial gauge A^ = 0. The main difference is as follows. In 
LFQCD, the finite energy density for physical states results explicitly in the asymptotic 
equation for transverse gauge fields in longitudinal infinity [see eq.(27)] due to the fact that 
the conjugate momenta of A l a are dependent variables in light-front form, [E l a = ^d + A l a 
is not a light-front time derivative). In axial gauge, A l a ,i = 1,2 and their conjugate mo- 
menta are all the dynamically independent variables. Thus, a similar condition proposed by 
Chodos leads to a very complicated formalism which may not be practically useful even for 
perturbation theory, as noted by himself |^2|. The second major difference is the vacuum. 
In axial gauge, the QCD vacuum is still complicated regardless of the boundary condition 
chosen. In such a case, it is very difficult (if not impossible) to do non-perturbative cal- 
culations before knowing the vacuum structure. In LFQCD, with antisymmetric boundary 
conditions, the vacuum is trivial and the nontrivial behavior of QCD would be manifested 
directly in Hamiltonian operators induced by the boundary integrals. Thus it is straight- 
forward to use quantum mechanical non-perturbative approaches to compute bound states. 
Moreover, in axial gauge, the boost invariance is not manifested kinematically so that it is 
not a good framework to study low-energy QCD, which deals with composite particles of 
quarks and gluons. In LFQCD, as we have mentioned in the introduction, boost invariance 
is a kinematical symmetry which is very useful in addressing hadronic structures. 

In summary, the essential point in determining nontrivial behavior of LFQCD seems to 
be the boundary conditions. A suitable choice of boundary condition for physical fields in 
LFQCD is crucial because it determines whether the nontrivial behavior of QCD can be 
decoupled from the vacuum so that the property of the trivial vacuum in LFFT becomes 
useful for solving hadrons from QCD. We have derived the canonical formulation of LFQCD 
with great care for boundary integrals, which have not been paid enough attention in 
previous investigations. We show that the boundary integrals are the source of the light- 
front linear infrared singularity and involve color charge densities and non-local behavior in 
the transverse direction that lead to non-local forces generated by the boundary integrals 
which are also responsible for the non-vanishing topological winding number. Clearly, our 
understanding of the physics from the boundary integrals in LFQCD is far from complete 
and much work remains to be done. Particularly, two questions are very interesting for 
the understanding of hadronic physics. One is which terms among the numerous non-local 
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interactions are essential for hadronic bound states. The other is how we can find an explicit 
field configuration for the non-zero winding number in LFQCD that satisfies the asymptotic 
behavior of eq.(27). These are two of the main problems in nonperturbative LFQCD for 
the future. 
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Appendix A. Canonical quantization of LFQCD 

A self-consistent formulation of LFQCD requires that the resulting Hamiltonian must gen- 
erate the correct equations of motion for the physical degrees of freedom (A i a ,iJ)+,ifi+). 
This appendix is devoted to the derivation of canonical quantization and a check of the 
consistency. 

To see how to correctly reproduce the Lagrangian equations of motion, we need to find 
consistent commutators for physical field variables. 

In the light-front gauge, the Lagrangian of eq.(7) is reduced to 

c = \d + Aid-Ai + - d~^U + ) - n. (A.i) 

The canonical momenta of the physical field variables A l a ,ip + ,ip+ are 

£ l a = \d + A\, ^ + = l -^ +) ^t_ = -iv+- ( A - 2 ) 

However, eq.(A.2) shows that all the canonical momenta are functions of the independent 
field variables. Thus, after determining all the Lagrange multipliers, the system is still a 
constrained Hamiltonian system. Usually, in order to quantize such a constrained Hamil- 
tonian system, one has to use the Dirac procedure, by imposing the so-called primary, 
second-class constraints E % a -\- ^d + A l a = (similarly for 7r^ + ^ ) to construct Dirac brackets. 
However, for these trivial primary constraints, the mathematically well-defined canonical 
one- form offers a rigorous phase space structure for canonical quantization [13|.0 In this 
appendix, we will use such an approach for light-front quantization. 



4 If the first-class constraints cannot be solved explicitly, for example for other gauge choices in light-front 
form, using the Dirac procedure may be necessary to construct the canonical commutation relations for all 
variables. 
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The phase space structure for the physical variables (£ \ , A l a ; ir^ + , ir^\ , ip + ) is deter- 
mined rigorously by rewriting eq.(A.l) as a Lagrangian one-form Cdx + (apart from a total 
light- front time derivative), 



Cdx + = ^2(£ l a dAi + Ti^ + d^ + + #J.?r^ - A ? a cZ£* - - V'+^O - ttdx + 

1 



^"Tafsdq 13 - Hdx + (A.3) 

where the first term in the right-hand side is called the canonical one-form of the physical 
phase space, and quark fields are anticommuting c- numbers (Grassmann variables). Cor- 
respondingly, the symplectic structure or the Poisson brackets of the phase space is given 
by 

u = ±r a pdq a dqf 3 or [<f ,<f\ p = T^. (A.4) 

Canonical quantization is realized by replacing the Poisson brackets by the equal-x + com- 
mutation relations 



[^,«°] = ir-J. (A.5) 



Explicitly 



1 



[Al(x),£l(y)} x+=y+ = i-5 ab 8»5\x - y), (A.6) 



.A, 



{ip+(x),Tr f+ (y)} x+=y + = i— S (x-y), (A.7) 
Hl(x),ir^ + (y)} x+=y+ = -i^-5 3 (x - y) (A.8) 



or 



[Al(x),d + Ai(y)] x+=y+ = i5 ab ^5 3 (x - y), (A.9) 
[Al(x),Ai(y)] x+=y+ = -iS ab ^]e(x- -y-)5 2 (x L -y ± ), (A.10) 



4 



{^ + {x)^\{y)} x+=y+ = A + 5-\x - y) (A.ll) 

where 5 3 (x — y) = S(x~ — y~)5 2 (x±_ — y±). All other commutators between the phys- 
ical degrees of freedom vanish. Note that, unlike in the instant form, the commutator 
[A l a (x), A 3 a (y)] does not vanish since the canonical momentum is a function of the coor- 
dinates. Eq.(A.lO) is defined consistently with the antisymmetric boundary condition, 
A l a (— oo) = — A l a (oc), where it is also required that 

lim e(x~-y~) = (A.12) 

Topologically, this requirement is identical to e(0) = 0. This can easily be understood if we 
divide the longitudinal line into boxes and define the e-function in each box. 
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Using the above basic commutation relations, it is straightforward to verify that the 
equations of motion are consistent with eqs.(2) and (3), 

= j[V>+ , H] 

= UgA~ - • (id± + gA±) + (3m} 

x / dx'~e(x~ - x'~){a± ■ {id± + gA ± ) + /?m} j (A.13) 

J — oo J 

1 Z" 00 ... 
= - / dx'-e(x--x , -)[Di b F^(x+,x'-,x 1 _))-D^(x+,x'-,x 1 _) 

- gf a (x + ,x'-,x ± ) - gr bc (A^d + Al)(x + ,x'-,x ± )} (A.14) 
where -D* b = ^afcc'* — gf abc A l c , and A~ and £7" are given by eqs.(19) and (15). 



Appendix B. Cancellation of linear infrared divergence 

In this appendix, we shall use the x + -ordering perturbative rule which we developed recently po| 
for the two-component LFQCD to check the cancellation of linear infrared divergence in per- 
turbative LFQCD. For the reader's convenience, we list some relevant diagrammatic rules 
in Table 1 for the following calculation. For a complete list of the x + -ordering perturbative 
rules and Feynman rules for two-component LFQCD, see Ref.[|f|. 



B-l. Tree level (qq scattering) 

The lowest-order qq scattering amplitude is given by 



M 



fi 



fi 



7*' 



(B.l) 



which corresponds to the diagrams shown in Fig.l. Using the rules listed above and the 
x + -ordering perturbative theory, we immediately obtain that 



M% + M° fi = g 2 T^ 3X lT q (p2,Pi)xixir q (P4,P3)X3 
f 1 1 



A-rl/ 



+ 



Pi -Pi -Pa + k Pi ~P2 -P3 

V^TJ^xi^xlxa NB, 
l 



2g 2 T$ 1 T&xtxiTmzxlx3 



WB 



where 



k l 



[Ar 1 



a-p2±-im i iCr-pijL+im 
— a — a 



h4 



[pi 



(B.2) 
(B.3) 

(B.4) 
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p i is the total energy of the initial state, = (pf — p2,p\ — p\i ^ Pl± + 2J -\ +m ). NB 

Pi ~P2 

denotes no boundary term and WB means including the boundary term. It follows that 
in the principal value prescription, MJ i + M« is free of infrared divergences, while 
without the boundary term has a 1/e 2 divergence when k + — > 0. When the boundary 
term is included, the 1/e 2 is cancelled [see (B.3)]. Therefore, it is necessary to include the 
boundary term in order to obtain a finite amplitude for the lowest-order qq scattering. A 



similar discussion for e + e scattering in LFQED is given in ref . 1 23 1 . 



B-2. Loop corrections (one-loop quark self-energy) 

Based on the x + -ordering perturbative theory, the quark on-shell self-energy (mass correc- 
tion) up to one-loop is determined by 

£(p 2 = m 2 ) = S a + S fc + S c . (B.5) 

The three terms in the right-hand side are denoted by the three diagrams shown in Fig. 2. 
Again, using the rules listed above, we find that 

2 n f dk + d 2 k ± 9(k + )6{p + - k+) j k l a-pj_-im i 

Sa " 5 Cf J [jfe+] \ 2 W\ jp+j ° 

i a ■ p± — a ■ k± + im 1 J k l a ■ p± — a ■ k± — im i 
~° [p+ - fe+i J [ |jfe+] [p+ - k+] ° 

_ ai cT-p ± + im } 1 



J p — A; — (p — k) 

S " " 25 W 16^ [A; + ]b+-fc+] (B - 7) 

E J 252c / J" ^Sg^ j kp+-W] - [( P ++k+r\ } NB (B8) 

\ 2 9 2 C f J dk ^J 1 { [p +\ +] i ~ [p ++Vp } WB 
where C/ = 4/3. A direct calculation shows that 

= ld?kJ\n4A (B.10) 



87T 2 p+ 7 VP 
#^/d 2 £; ± (-2) NB 

which tells us that in the one-loop correction to the quark self-energy, the one-gluon 
exchange contains both linear and logarithmic infrared divergences. The instantaneous 
fermion interaction contains only one logarithmic divergence (see in eq.(B.lO)), which 
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cancels the logarithmic diverg ence in S a . The naive instantaneous- gluon interaction (namely 
the linear potential in the longitudinal direction) is free of infrared divergence. Therefore, 
without boundary term, the quark mass correction involves a linear infrared divergence, 
which is an inconsistent solution. By combining the boundary term with the linear po- 
tential, we see that S c has a linear infrared divergence which precisely cancels the same 
divergence in S a . Thus the quark mass correction is now free of infrared divergences, 

2 A 2 

5m 2 =p+S = ^rCf In — = + finite (B.12) 

where A is the transverse momentum cut-off. In eq.(B.9), the coefficient (1/4) in the mass 
correction is different from the covariant result (3/8) because the regularization scheme is 
different. This coefficient is the same as that in the light-front calculation with dimensional 
regularization in the transverse direction and the explicit cutoff in the longitudinal direction 



[41]. Note that in their calculation, the expressions of eqs.(B.9-ll) are different but the 



sum is the same as eq.(B.12) where the linear divergence is also cancelled. 
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Table 1: Some x + -ordering diagrammatic rules for the two-component LFQCD: 



a 



Pi 




P 



P2 



-9K 1 )Xj3 \ JFFJ uTi — CT _CJ — r^Fi — f X^ 



Pi 



a . 



fci 



j juuui 



A: 2 



P2 







Pi P2 a 

a t » p 



7 



P3 



PA 



2g 2 (T a T a )xl Xa 

2g 2 (T a T a )x i f3Xa ^ 



Pi. 



NB 
WB 



Remarks. The Xa are two-component light-front quark spinors, XT = (o) anci Xl = (?)> 
and the e % are two-component gluon polarization vectors. We have also used the Majorana 
representation of the 7 matrices in the above realization of the two-component formulation. 
In additional, the internal integral for quark momentum is given by / dp 1 ^ ± 6{p + ), while 

for gluon it is / ^j^r^^ry- The light- front free quark and gluon energies are p~ = P± ^™ 
and k~ = 
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Fig.l. The x + -ordering graphs for the lowest-order qq scattering in perturbative LFQCD. 
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Fig. 2. The x + -ordering graphs for the one-loop correction of quark self-energy in pertur- 
bative LFQCD. 
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